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Abstract 



An isospin analysis of B ^ tttt decays yields sin 2a, where a = 
arg[— and Vij is a CKM matrix element, without hadronic 
uncertainty if isospin is a perfect symmetry. Isospin, however, is broken not 
only by electroweak effects but also by the u and d quark mass difference. 
The latter drives vr*^ — r], rj mixing and converts the isospin-perfect triangle 
relation between the B tttt amplitudes to a quadrilateral. The combined 
isospin- violating effects impact the extracted value of sin 2a in a significant 
manner, particularly if the latter is small. 
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In the standard model, CP violation is characterized by a single phase in the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix, rendering its elements complex. Although CP violation 
has been known in the neutral kaon system since 1964, the absence of definitive evidence 
for a non-zero e' parameter leaves the above standard model picture unsubstantiated [|l|. 
Indeed, probing the precise mechanism of CP violation will be the primary mission of the 
future B factories. The CKM matrix of the standard model is unitary, so that determining 
whether or not this is empirically so results in a non-trivial test of the standard model's 
veracity. In the CKM matrix, only one combination of rows and columns results in an 
unitarity test in which all the terms are of the same approximate magnitude |0 ; this is the 
unitarity triangle 0. Empirically determining whether its angles, termed a, (3, and 7, sum 
to vr and whether its angles are compatible with the measured lengths of its sides lie at the 
heart of these tests of the standard model. 

In the decay of a neutral B meson to a CP eigenstate /cp, CP violation can be generated 
through B^-^ mixing, specifically through the interference of B^ ^ fcp and B^ B° 
fcp- Thus, weak phase information can be extracted from the time-dependent asymmetry 
A[t), defined as 

^ nB\t) ^ fcp) - Tj^jt) ^ fcp) 
T{BO{t)^fcp) + T(B\t)^fcp)' 

noting B°{t = 0) = 5° and B {t = 0) = B . Indeed, were only amplitudes with a single 
CKM phase to contribute to = 0) ^ fcp, the weak phase information could be 

extracted directly from A{t) without hadronic ambiguity 0. Unfortunately, however, either 
penguin contributions or a plurality of tree-level contributions arise to cloud the above 
analysis 0]. 

Nevertheless, the quantity sin 2a, where a is the usual CKM angle a = 
aTg[—VtdV^l/ (yudV*f^)] 0, can be extracted without penguin "pollution" from an isospin 
analysis of B ^ titt decays if isospin is a perfect symmetry [Q. In this limit, the Bose sym- 
metry of the J = vrvr state permits amplitudes merely of isospin / = 0, 2. This implies that 
the amplitude B^ ti^tt^ is purely 1 = 2. Thus, as two independent amplitudes describe 
the three amplitudes B^ tt^tt , B tt^tt^, and B they can be drawn as a 

triangle. A triangle can also be formed from the amplitudes B^ tt^tt^, B tt+tt^, and 
B^ TT^ir^, with the B^ vr^vr" amplitudes forming a common base. The strong penguin 
contributions are of A J = 1/2 character, so that they cannot contribute to the 1 = 2 am- 
plitude and no CP violation is possible in the vr^vr'' final states. This implies that the CP 
violation due to the penguin contribution in B^ tt+tt", or analogously in i? tt+tt^, 
can be isolated and removed by identifying the relative magnitude and phase of the 1 = 
and 1 = 2 amplitudes. 

It is our purpose to examine the manner in which isospin-violating effects impact the 
extraction of sin 2a as determined in i? — >■ tttt decays 0. In the standard model, isospin 
is an approximate symmetry. Isospin is broken not only by electroweak effects but also by 
the strong interaction through the u and d quark mass difference. Both sources of isospin 
violation generate A J = 3/2 penguin contributions, but the latter also drives vr'' — 77,77' 
mixing 0, admitting an / = 1 amplitude [0. These latter contributions convert the triangle 
relations between the amplitudes to quadrilaterals. The effect of electroweak penguins has 
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been studied earlier in the literature, and is estimated to be small [l^,^,^]. Nevertheless, 
when all the effects of isospin violation are included, the impact on the extracted value of 
sin 2a is significant. 

To review the isospin analysis in i? — > tttt decays 0], let us consider the time-dependent 
asymmetry A{t) 0: 

where r;„ = (V,tV„/V,,V,;)(A''^- /A'''-) = e-"*-^ and yt'" = A(B« ^ fcp). We 
assume the mass eigenstates Bl and Bh have the same width and a mass difference Am = 
Bh — Bl. The sin(Amt) term, resulting from B'^-B mixing, is linear in rj^p and thus 
is of especial interest. If fcp = vrvr, then the presence of penguin contributions implies 
^fcp ^ -jfcp^ denote the amplitudes B+ tt+ti^, B° 7r°7r°, and 5° Tr+Tr", by 
y4+°, and respectively, and, following Ref. [Q, we write 

^A+-=A2-Ao ■,A'' = 2A2 + Ao ; = 3^2 , (3) 
noting analogous relations for A~^, A^^ , and A'^ in terms of A2 and Aq. Thus, 

(1-.)' 

where z(z) = ^0/^2(^0/^2) and A2/A2 = exp(-2i0t) with 0f = arg(KdK*6) and (pm + (pt = 
/3 + 7 = 7r — ain the standard model [Q. Given l^"*""], |^°°|, l^"'"''!, and their charge 
conjugates, the measurement of Imr7r+7r- determines sin 2a, modulo discrete ambiguities in 
arg((l — z)/(l — 2;)). The latter can be removed via a measurement of ImrTrO^o as well [Q. 

We proceed by computing the individual amplitudes using the AB = 1 effective Hamil- 
tonian resulting from the operator product expansion in QCD in next-to-leading logarithmic 
(NLL) order [0,0], using the factorization approximation for the hadronic matrix elements. 
The factorization approximation, which assumes the four-quark-operator matrix elements 
to be saturated by vacuum intermediate states, finds theoretical justification in the large Nc 
limit of QCD [|ri|] and phenomenological justification in comparison with empirical branching 
ratios [0; nevertheless, it is heuristic. We adopt it in order to construct concrete estimates 
of the effects of isospin violation in the decays of interest. In this context, we can then 
apply the isospin analysis delineated above to infer sin 2a and thus estimate its theoretical 
systematic error, incurred through the neglect of isospin violating effects. 

The effective Hamiltonian for b ^ dqq can be parametrized as [|l^ 



V2 



V^,V:,iC^O^ + C^Ol^) + Ki,V:,{C,Ot + C^O^^) - VtiMt^d E C^0^ + CgO, 



(5) 



where Oi and Og are as per Ref. |]T0|; we also adopt their Wilson coefficients Ci and C^, com 



puted in the naive dimensional regularization scheme at a renormalization scale of /i = 2.5 
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GeV [0]. In NLL order, the Wilson coefficients are scheme- dependent; yet, after comput- 
ing the hadronic matrix elements to one-loop-order, the matrix elements of the effective 
Hamiltonian are still scheme- independent |]13|. This can be explicitly realized through the 
replacement {dqq\n'''\b) = (Gp/v^)(f/gg|[K6Krf(Q''0]^+QffO^)-l^tb\/;^E-£3Q''C',]|6)*'--, 
where "tree" denotes a tree- level matrix element and the Cf^ are from Ref. iTOll. The Cf^ 



are complex and depend on both the CKM matrix parameters and fc^, where k is the 



momentum transferred to the qq pair in 6 — »• dqq. Noting Ref. 0] we use p = 0.12, rj = 0.34, 
and A = 0.2205 |lT0| , p!5| unless otherwise stated. One expects m^/4<fc^<m^/2 |jl6| 



we use 



k'^/ml = 0.3, 0.5 in what follows. 

To include the effects of tc^-i], rj' mixing, we write the pion mass eigenstate Ivr") in terms 
of the SU{3)f perfect states |03) = \uu — dd) / \/2, l^g) = \uu + dd — 2ss) / \/6, and |0o) = 
\uu + dd+ ss) I \f?). To leading order in isospin violation |^ 

k°) = 103) + £(cos^|08) - sin^|0o)) + £'(sin^|08) + cos^|0o)) , (6) 

where |?7) = cos6'|08) — sin6'|0o) +0(£), and \r]') = sin6'|08) + cos6'|0o) +0{e'). Expanding 
QCD to leading order in l/N^, momenta, and quark masses permits the construction of a 
low-energy, effective Lagrangian in which the pseudoscalar meson octet and singlet states 
are treated on the same footing Diagonalizing its quadratic terms in 03, 08, and 0o 

determines the mass eigenstates 7r°, rj, and r]' and yields e = eoxcosO and e' = — 2eoXsin^, 
where x= l+(4m^-3m^-m^)/(m^-m^) ^ 1.23, x = l/X;£^o = V^imd-mu)/{4:{ms-m)), 
and 771 = {mu+md) /2 0. Thus the magnitude of isospin breaking is controlled by the SU(3)- 
breaking parameter ms — rh. The rf-r]' mixing angle 6 is found to be sin 2^ = —{A\/2/2>){rri]^ — 
m^)/(m^, — m^) ^ —22° 0. The resulting e = 1.14eo is comparable to the one-loop-order 
chiral perturbation theory result ofe = 1.2360 in ^ tt+tt^tt^ [T^JsI]. Using mq{fi = 2.5 GeV) 
from Ref. [ig], we find e = 1.4 ■ lO'^ and e' = 7.7 ■ IQ-^. 

We now compute the matrix elements of the above effective Hamiltonian in the cases of 
interest. We define the decay constants (7r~(p)|(i7^75M|0) = — z/,rP/i and (0j(p)|M7^75tt|O) = 
—if^jPfj,, and use SU(3)/ to relate to Finally, using the quark equations of motion 
with PCAC and introducing = Cf + Q^JN^ for i odd and = Cf + CfjN^ for i 
even, the B~ 7r~03 matrix element in the factorization approximation with use of the 
Fierz relations [01 is 



{niu + md){mb ~ iTLu) 

,3 1 m2„(a6 - \as) 

X (a4 + -(a, - a,) - -a,, + ^^^^^ _ ))] • 

The transition form factors are given by FB-Tviq"^) = {^"b- " '"^^-)-^(f^'^(0)/(l ~~ 
where we use F^^^{0) = 0.33 and Mo+ = 5.73 GeV as per Refs. Also ^5^3 = 

Fbtt/V^, Fb(I)s = Fbtt/VQ, and -Fb^q = Fb-k/V^- In the presence of isospin violation, the 
7r^03 amplitude is no longer purely 1 = 2. However, if vr^-r/, rj' mixing is neglected, 
A^' + 2A = V2A , from Eq. |, is still satisfied, as we ignore the small mass differences 
m^± — m^o and m^i — m^o. Now with n^-rj, rj' mixing. 
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= (7r-03|W''^|5-) + £8(vr-08|^'^|5-) + £o(vr-0o|^''^|5-) (8a) 

= (0303|^^''|5") + 2£8(0308|^'=^|5°) + 2^0 (0300 1^°) , (8b) 



where eg = £cos9 + e' sm6, Eq = e'cos6' — £sin6', and further details appear in Ref. |]T9[. 
The B ^ TiTT amphtudes now satisfy 

~A^' + 2A°° - V2 = 458(03081^''' 1 5°) + 4^0 (0300 1 ^''^ 1 5^) 

- V2£8(vr-08|H^^|5-) - V2£o(vr-0o|H'''|5-) , (9) 

and thus the previous triangle relation becomes a quadrilateral. Numerical re- 
sults in the factorization approximation for the reduced amplitudes and A^, 

where ^ 2A°7((G^/v^)^K6Kd), = A+7((G^/v^)^K6Kd), and A],^ ^ 

^/2A-°/{{GF/V2)iVubV*a), with = 2, 3, cx) and P/m^ = 0.3,0.5 are shown in Fig. |l|. 
Here the parameter Nc defines different hadronic models; iVc = 2, 3 bound the value favored 
from fits to measured branching ratios [0. A'j^^ and A]i^ are broken into tree and penguin 



contributions, so that y4j° = T^+(/,3 + Ptt+tt" and A]^^ = T^,-^^ + P-n-T^o, where PjriTrO is de- 
fined to include the isospin-violating tree contribution in A^^ as well. The shortest side in 
each polygon is the vector defined by the RHS of Eq. For reference, note that the ratio 
of penguin to tree amplitudes in n-n^ is |P|/|T| ~ (2.2 - 2.7)%\VtbV;a\/\Vy,bV*a\ for 

A^c = 2, 3 and k"^ above. Were electroweak penguins the only source of isospin violation, then 
\P\/\T\ ~ {lA-1.5)%\VtbV;a\/\VubV*a\, commensurate with the estimate of 1.6% in Ref. [fl. 
Following Ref. 0], the determination of Imr^+^- yields sin 2a, modulo discrete ambiguities 
in z and z as in Eq. ^. However, sin 2a can be determined uniquely through a comparison 
with sin 2a from Imrj^o^o, as only one pair of the sin 2a extracted from Imr^+^~ and Imr^o^o 
likely match. As |A2|/|v42| ^ 1 in the presence of isospin violation, we have retained this 
explicit factor in the last equality of Eq. ^ in order to extract sin 2a as accurately as possi- 
ble. The values of sin 2a extracted from the amplitudes in the factorization approximation 
with Nc and k'^/ml = 0.5 are shown in Table | for a variety of input values of sin 2a — 
the results for k'^/ml = 0.3 are similar and have been omitted. In the presence of 7i^-ri,7]' 
mixing, the A^ , A^", and A^ amplitudes obey a quadrilateral relation as per Eq. so 
that the amplitudes of interest need no longer form triangles. Consequently, the values of 
sin 2a extracted from Imr^+^- and Imr^o^o can not only differ markedly from the value of 
sin 2a input but also need not match. The incurred error in sin 2a increases as the value 
to be extracted decreases; the structure of Eq. || suggests this, for as sin 2a decreases the 
quantity Im((l — z)/{l — z)) becomes more important to determining the extracted value. 
It is useful to constrast the impact of the various isospin-violating effects. The presence of 
AI = 3/2 penguin contributions, be they from m„ ^ rrid or electroweak effects, shift the 
extracted value of sin 2a from its input value, yet the "matching" of the sin 2a values in 
Imr^+^~ and ImrT^OT^o is unaffected. The mismatch troubles seen in Table | are driven by 
TT^-f], Vj mixing, though the latter shifts the values of sin 2a in \m.r.^+T^~ as well. Picking the 
closest matching values of sin 2a in the two final states also picks the solutions closest to 
the input value; the exceptions are noted in Table |. The matching procedure can also yield 
the wrong strong phase; in case b) of Table | with = 3, oo, the triangles of the chosen 
solutions "point" in the same direction, whereas they actually point oppositely. 

If l^^^l and are small the complete isospin analysis may not be possible, so that 
we examine the utility of the bounds proposed in Ref. [EO] on the strong phase 25true = 
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arg((l — z)/{l — z)) in Eq. ^. The bounds 2(5gqi and 2(5gqii given by Eqs. 2.12 and 2.15, 
respectively, in Ref. follow from Eq. and thus can be broken in the presence of isospin 
violation. As shown in Table |, the bounds typically are broken, and their efficacy does not 
improve as the value of sin 2a to be reconstructed grows large. 

To conclude, we have considered the role of isospin violation m B ^ -k-k decays and have 
found the effects to be significant. Most particularly, the utility of the isospin analysis in 
determining sin 2a strongly depends on the value to be reconstructed. The error in sin2Q; 
from a ImrTr+jr- measurement can be 50% or more for the small values of sin2Q; currently 
favored by phenomenology [^,^,^; however, if sin 2a were near unity, the error would 
decrease to less than 10%. The effects found arise in part because the penguin contribution 
in TT+TT", e.g., is itself small; we estimate |P|/|T| < 9%|\4bV^^|/|Kf)Kdl- Relative to 

this scale, the impact of tt^-t], t]' mixing is significant. Yet, were the penguin contributions in 
B ^ irn larger, the isospin-violating effects considered would still be germane, for not only 
would the AI = 3/2 penguin contributions likely be larger but the B ^ ttt] and B nt]' 
contributions could also be larger as well ||22|. To conclude, we have shown that the presence 
of TT^-i], 1]' mixing breaks the triangle relationship, Eq. ^ usually assumed and can mask 
the true value of sin 2a. 

S.G. thanks S. Cotanch for a suggestive query and D. London, T. Cohen, A. Kagan, 
A. W. Thomas, and M. Wise for helpful discussions. This work was supported by the DOE 
under DE-FG02-96ER40989. 
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TABLES 



TABLE L Strong phases and inferred values of sin 2a from amplitudes in the factorization 
approximation with Nc and k"^ /ml = 0.5. The strong phase 2(5true is the opening angle between 
the Ar^ and Ar' amplitudes in Fig. 1, whereas 2(5gl is the strong phase associated with 
the closest matching sin 2a values, denoted (sin2a)GL! from Imr^+^-/Imr^o^o, respectively. The 
bounds |25gqi| and |2(5GQn| on 25true from Eqs. 2.12 and 2.15 of Ref. ||2^ are also shown. All angles 
are in degrees. We input a) sin 2a = 0.0432 0,|l|, b) sin 2a = -0.233 (p = 0.2, r? = 0.35) |2l|, 
and c) sin 2a = 0.959 [p = —0.12). * The matching procedure fails to choose a sin 2a which is 
as close to the input value as possible. ^ The discrete ambiguity in the strong phase is resolved 
wrongly. 



case 






|2<5gqi| 


|2(5gqii| 


|25gl| 


(sin2a)GL 


a 


2 


24.4 


26.1 


15.8 


16.6 


-0.0900/-0.0221 


a 


3 


24.2 


16.9 


16.1 


16.2 


-0.0926/0.107 


a 


oo 


23.8 


59.4 


25.1 


23.6 


0.0451/0.394 


b 


2 


19.6 


23.4 


12.1 


12.9 


-0.343/-0.251 


b 


3 


19.4 


13.5 


12.9 


13.0 


-0.719/-0.855(*) 


b 


CO 


19.2 


59.9 


23.6 


0.76 


-0.550/-0.814(*,t) 


c 


2 


28.3 


36.5 


20.4 


21.0 


0.917/0.915 


c 


3 


28.0 


24.0 


19.1 


19.0 


0.905/0.952 


c 


oo 


28.3 


36.5 


20.4 


21.0 


0.917/0.915 
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FIGURES 

0.50 I ■ 1 ■ 1 




Re{A,} (GeV') 



FIG. 1. Reduced amplitudes in S — tttt in the factorization approximation with [Nc, k^/ml] 
for a) [2,0.5], b) [3, 0.5] (sohd hne) and [3, 0.3] (dashed hne), and c) [oo,0.5]. 
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